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We investigate numerically the statistical properties of spectra of two-dimensional 
disordered systems by using the decimation method applied to the Anderson model. 
Statistics of the spacings calculated for system size up to 1024 x 1024 lattice sites 
exhibits a crossover between Wigner and Poisson distributions. We perform a self- 
contained finite-size scaling analysis to find a single-valued one-parameter function 
7(L/£) which governs the crossover. The scaling parameter £,(W) is deduced and 
compared with the localization length. j(L/£) does not show critical behavior and 
has two asymptotic regimes corresponding to weakly and strongly localized states. 
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The single-parameter scaling theory of Anderson localization suggests that there exists 
no conducting state in an infinite two-dimensional (2D) disordered system at zero tem- 
perature [Q], since all one-electron wave functions are localized even for arbitrarily small 
fluctuations of a random potential @. The corresponding energy levels are distributed com- 
pletely randomly. One expects that there are no spectral correlations at any energy scale. 
Therefore, for any non- vanishing disorder the distribution of neighboring level spacings P{s) 
is given by the Poisson law Pp(s) = exp(— s) in the thermodynamic limit. (Energy is mea- 
sured in units of the mean level spacing A.) 

For finite systems P(s) ~ Pp{s) only if the disorder W is strong enough such that the 
localization length £ is shorter than the system size L. Indeed, for fixed W^/fl the level 
spacing distribution was numerically shown to approach Pp{s) with increasing L 0. With 
decreasing W the electron states change gradually from "strongly localized" to "weakly 
localized", and the energy levels become more correlated. When £ ^> L, their statistics is 
well described by the Wigner-Dyson theory for eigenvalues of random matrices [Q . 

As a consequence of the transition from strong to weak localization, P(s) is expected 
to undergo a continuous crossover between Pp(s), and the Wigner surmise Pw{s) = 
txs/2 exp(— 7rs 2 /4). In the weakly localized and non-diffusive regimes the spectral fluctu- 
ations were studied in terms of the two-level correlation function by the use of sophisticated 
analytical techniques Small corrections to the Wigner-Dyson statistics were obtained 

in 0. Although the absence of the metal- insulator transition for statistics of the spectra 
in 2D-system was already mentioned by Shklovskii et al ||, so far extensive investigations 
in the intermediate regime, where the deviations from both of the limiting distributions 
are considerable, have not yet been carried out. An important question is which scaling 
law governs the crossover of the level statistics. For three-dimensional (3D) disordered sys- 
tems such a scaling law shows critical behavior at a finite disorder which is the signature 



of a metal-insulator transition [|3|-|T0|j. If no derealization transition exists in 2D [TT[, P(s) 
should not exhibit critical behavior at any finite disorder. In this paper we provide strong 
numerical evidence that this is indeed the case. 



The scaling behavior of P(s) with L and W is determined by using the results for the 
energy spectrum that were obtained by a combination of the decimation method and a 
Lanczos procedure. Due to the achieved large system sizes, the scaling analysis could be 
performed very accurately Independently of other methods, we calculate a one-parameter 
scaling variable j(L/£) related to P(s) which governs a crossover between the two above 
limits. It turns out to be monotonic without any critical point W c , for which j(L, W c ) = 
const, that would indicate a derealization transition. A similar variable ^{L/^tmm) was 
claimed in || to behave as a single- valued function, where £tmm is the localization length 
obtained from the transfer- matrix method (TMM) ||12|| . We find that in the strongly localized 
regime the deviation from the Poissonian is described by a linear law 7 oc L/£. In the weakly 
localized regime 7 tends towards the Wigner limit. From spectral statistics we determine 
the scaling parameter £, which can be interpreted as the localization length. Its dependence 
on W is consistent with the results of previous calculations using the TMM [0,0 • 

We consider the Anderson model (AM) with diagonal disorder H = J2i e i I O ( i I 
+tJ2<i,j> I i ){ 3 I) where i labels the L 2 sites of a square lattice of the linear size L. 
The on-site energies £j are independently distributed at random within an interval of width 
W. The second sum is taken over all pairs of nearest-neighbor sites. It corresponds to 
the kinetic energy. For the calculations we used a decimation algorithm described in PjMl, 
which enables us to investigate the level statistics near the band center E for systems 
containing up to 10 5 -10 6 lattice sites. The algorithm consists of a sequence of iterative 



steps [[T4J], in which the spectrum of a system is computed by using the amplitudes of 
wave functions and corresponding eigenvalues of four smaller systems of half of the linear 
size. Each step consists of (i) a unitary transformation of the secular matrix to the basis 
set of the smaller lattices, (ii) a truncation of the corresponding Hilbert space to a fixed 
number m of states near E — 0, and (iii) a diagonalization of the truncated Hamiltonian. 
The iteration starts with the original 2D lattice of the AM. The latter is divided into many 
independent small square cells of linear size L D , which is defined by the disorder. Subsequent 
larger squares are constructed by coupling four smaller squares via hopping elements t only 



between corresponding boundary sites. In the last step of the iteration, periodic boundary- 
conditions are imposed. Using ensembles of up to k ~ 1000 different realizations of the 
random potential, we calculate P(s, L, W) at E . 

As expected, we find the Poisson distribution for large L and large W and the Wigner 
distribution for small W and small L. The crossover between the two limiting cases can be 
detected for fixed disorder with increasing L, or for fixed L with increasing W. Figure [l] 
shows that for W = 4 the level spacing distribution is only slightly different from the 
Wigner surmise. On the other hand, P(s) for W = 7 already resembles strongly the Poisson 
distribution. This can be explained by the considerably different values of the localization 
length for these two disorders. In the former case, £ ^> L while in the latter case £ <C L (see 
the inset of Fig. 3 below). In order to test our method we also calculated P(s) by using 
a Lanczos procedure fl~5| , especially optimized for large sparse matrices |10|. We were able 



to compute the exact spectrum for systems up to L 2 = 256 2 . The results of the decimation 
method for P(s) are in good agreement with those obtained by the exact diagonalization 
(inset of Fig. [I]). For larger systems the Lanczos procedure failed even for a very narrow 
energy interval due to limitations of computing time, while we could enlarge the system up 
to L 2 ~ 10 6 lattice sites by using the decimation algorithm. The striking advantage of the 
latter is that it allows to go beyond the limits of conventional diagonalization procedures, 
at the expense of being restricted to a small (and constant) number of energy levels around 
the energy of interest. 

The decimation scheme is applicable as long as V <C A, where V is the typical matrix 
element between states of the coupled blocks (see fljbj). Therefore one should start the 
calculations with such L Q for a given W, so that this condition is fulfilled. The choice of m 
defines the relative error of the retained eigenvalues SE/A ~ V 2 / (A 2 m) <C 1 resulting 
from the neglection of the discarded levels. The errors of the consecutive iterations steps are 
mostly dominated by those of the first step, which diminish with decreasing coupling, i.e. 
with increasing W. In order to control the precision of P(s) we choose L Q and m depending 
on W so that after several iterations V/A ^ 0.01. The initial size grows from L a = 16 for 



W > 8 to L Q = 256 for W < 5 and m ranges typically from 25 to 100. For the sizes smaller 
than L Q the exact diagonalization is applied. 

In order to investigate the scaling behavior of the level statistics quantitatively, we con- 
sider the integral I(s) = Jq P(s')ds'. It describes the probability to have a spacing smaller 
than s. Numerically, I(s,L,W) can be obtained with a higher accuracy than P(s,L, W). 
Figure @ shows how I(s) changes from Iw( s ) = 1 — exp(— its 2 /A) to Ip(s) = 1 — exp(— s) for 
fixed disorder in approaching the thermodynamic limit. For example, the small- s behavior 
of I(s, L, W) changes from quadratic for small to linear for larger system sizes. Also, the 
probability to find large spacings is enhanced, when increasing L. 

The next step is to determine the deviation of the level statistics from the Poissonian 
limit due to the finiteness of the system. We introduce the function 

7 (L, W) = ± £ [1(8, L, W) - I P (s)] ds, (1) 

which is similar to the quantities considered in Ref. |8]^T0|. The normalization constant iV 
is chosen such that 7 = 1, if J(s) = I\y(s). For the scaling analysis the choice of the upper 
limit So of the integral does not play an important role [|TUJ. By optimizing the numerical 
accuracy we fixed it to so = 1.25A. For an infinite 2D disordered system 7 = for all 
W 7^ 0. For finite systems, 7 is expected to change continuously from unity to zero with 
increasing L ||16|| . 

We computed 7 for W from 4 to 12 and L from 8 to 1024 (inset of Fig. 0). One observes 
that 7 decreases with increasing L for all W . For weak disorder it is close to unity and 
decreases slowly, as the system doubles in size, while for larger disorder the decrease is more 
rapid. Correspondingly, for a fixed L, 7 diminishes as the disorder becomes stronger. The 
data for W •> 7 saturate towards 7 = 0. Altogether, 7(L, W) undergoes a smooth monotonic 
crossover from unity to zero as L grows, with no critical behavior, at least for the disorders 
considered here. This is markedly different from the 3D case, where a critical disorder W c 7^ 
exists, below which the electron states are delocalized, and 7 should increase with L PJIUfl. 
We have estimated the error of 7, £7/7 ~ {$1 /l) s tat + {$l/l)syst- The former is of order of 
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{kml(s ))~ l l 2 < 2%. The latter can be roughly estimated as I(s )/(N r y) 5E/A and is about 
0.5% near L D /£ ^ 1. In the strong localization it decreases as 0.02 L D /£ exp(— L /£). Thus, 
the statistical error due to the finite number of samples is always larger than the systematic 
one resulting from the approximation. 

Assuming that 7 obeys a one-parameter scaling law j(L, W) = f(L/£), one can find 
the scaling parameter £,(W) such that all numerical data fall onto a common curve. We 
fixed the data for W = 7 and then shifted the data corresponding to other W to the left 
and to the right in a logarithmic scale, until the set of data of a given W overlapped for 
at least one value of L with those corresponding to another W. One can determine £ by 
this procedure up to a common numerical factor. The latter corresponds to a shift of all of 
the data in the inset of Fig. ||] by the same amount along the x-axis. We found this factor 
by adjusting £ at large W to the localization length, which was earlier found by using the 
TMM ||13|| . For this fit we choose ^tmm(^ = 10) which was smaller than the width of the 
quasi- ID strips used in the TMM calculations. As shown in Fig. |3], the data of 7(L/£) for 
all pairs of {L, W} collapse onto an one-branch curve within the numerical accuracy. Thus, 
our calculations confirm by a completely independent procedure that the crossover from 
the weak to the strong localization in 2D, which corresponds in the language of the level 
statistics to the crossover from the Wigner surmise to the Poissonian distribution, is governed 
by one-parameter scaling. The left, very flat part of the curve in Fig. |3], obtained for small 
disorders, corresponds to the weak localization. Here, the data for In 7 approach zero because 
the scaling parameter £,(W) becomes large compared to the size of the system. The small, 
but increasing deviation from 7 = 1 for L/£ -C 1 is in agreement with analytical results 
obtained earlier for the corrections to the Wigner-Dyson statistics in the case Sq ~ A 
Details of a quantitative comparison which involves also a calculation of the conductance of 



the system are given elsewhere |17 



When L/£ ^> 1 the data are well described by 7 = 7 (£/L), where 7 Q ~ 1.9. An 
analogous linear behavior was found previously for the Lyapunov exponents [|1^]. The fact 
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that for strong localization the conductance g oc exp(— 2L/£) suggests that 7 oc — (lng) -1 
for L/£ ^> 1. The intermediate part of the graph, where L ~ £, shows the smooth crossover 
between the two limits. Since 7(L/£) consists of only a decaying branch, we conclude that 
there is no metal-insulator transition at least for W > 4. The inset of Fig. shows £(W). 
For W < 6 one observes a faster increase of £ with decreasing than the localization length 
computed via TMM in jn|. On the other hand, we find good agreement with the values of 



Ctmm, reported recently in Ref. [0. It is reasonable to idendify our scaling parameter £ with 
the localization length of the electron states. The disorder dependence of £ is consistent with 



the theoretical prediction of an essential singularity of the localization length at W — [18 



in the form £(W) = aexp[(W /W) 2 }, where the length scale a is of order of the lattice 
constant. By fitting the calculated data of £ to this formula we found a = 1.1 ± 0.25 and 
Wo = 11 ± 3. 

In conclusion, we performed a finite size scaling analysis of the fluctuations in one- 
electron spectra of 2D disordered systems by using a combination of exact diagonalization 
and decimation algorithm. It enables us to treat systems up to 1024 2 lattice sites. The 
nearest-neighbor level spacing statistics described by the variable 7 (Eq. [I]) indicates a 
smooth crossover from strong to weak localization when the disorder and the system size 
decrease. The crossover is found to obey a one-parameter scaling law which does not show 
critical behavior. This confirms the earlier conjecture |l[] that in 2D there are no true metallic 
states in the thermodynamic limit. Our present results for the scaling parameter £ extracted 
from the level statistics are quantitatively consistent with those obtained previously for the 
localization length by using a completely different numerical approach |13| , and conclusions 
based on perturbative calculations |Tj| . 
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FIGURES 

FIG. 1. The probability density of neighboring spacings P(s) for system size I? = 256 2 at 
different disorders. Solid and dotted lines: Wigner surmise and Poisson law, respectively. Inset: 
P(s) at W = 6 obtained by the decimation and the Lanczos method. 

FIG. 2. Total probability distribution of spacings I(s) at the disorder VF=6 for various 
sizes. Solid and dotted lines: Iw( s ) an d -^p( s )> respectively. Inset: j(L) for different disorders 
W = 4.0 (e), 4.5 (★), 5.0 (o) , 5.5 (+), 6.0 (•), 7.0 (A), 8.0 (x), 9.0 (□), 10.0 (y), 12.0 (o). 

FIG. 3. The parameter 7 as a function of L/£(W). Straight line: limit of strong localization. 
Symbols for W are the same as in the inset of Fig. 0. Inset: disorder dependence of the scaling 



parameter £ (•) along with the results obtained for the localization length by TMM |12| (x) and [13] 
(A). 
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